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INTRODUCTION 


We  are  investigating  a  new  approach  to  detection  of  breast  cancer  using  optical  vascular 
function  imaging.  Research  on  cancer  therapy  has  revealed  unusual  properties  of  tumor 
vasculature  produced  through  angiogenesis  that  have  high  promise  for  breast  cancer  detection. 
We  are  using  this  unusual  behavior  as  a  contrast  mechanism  by  performing  optical  imaging 
during  administration  of  different  levels  of  oxygen  and  carbon  dioxide.  We  have  constructed 
imaging  systems  using  immersion  in  a  tissue  phantom,  improved  the  sensitivity  and  speed  of  the 
imaging  systems,  used  mass  flow  controllers  to  regulate  inhalation  gas  compositions,  imaged  5 
different  cancer  models,  and  compared  dynamic  and  static  measurements.  Our  results  support 
our  hypothesis  that  the  additional  contrast  from  imaging  vascular  function  is  superior  to  static 
(conventional)  optical  imaging.  The  results  of  this  IDEA  project  have  been  the  basis  of  new 
funding  for  clinical  testing.  We  believe  that  this  approach  can  form  the  basis  for  a  cancer 
imaging  modality  combining  high  sensitivity  and  specificity.  New,  sensitive,  specific,  and  low 
cost  imaging  technologies  are  very  important  for  detection  and  treatment  of  breast  cancer,  and 
should  lead  to  reduced  morbidity  and  mortality  from  this  disease. 


BODY 

The  following  sections  will  briefly  summarize  the  accomplishments  on  each  of  the 
project  tasks. 
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TASK  1.  MODIFY  AND  CONSTRUCT  IMAGING  SYSTEMS 


Two  optical  systems  were  used  in  our  studies.  These  include  a  frequency  domain  system 
for  quantification  of  scattering  and  absorption,  and  a  CW  system  for  monitoring  rapid  changes  at 
different  wavelengths.  Work  performed  on  these  systems  is  summarized  below. 

The  frequency  domain  apparatus  developed  in  our  laboratory  measures  the  optical  delay 
of  light  through  a  medium  in  order  to  quantitatively  determine  absorption  and  scattering 
properties  on  human  or  animal  tissue.  By  measuring  the  phase  shift  and  amplitude  decay  of  the 
modulated  light,  one  can  extract  the  scattering  and  absorption  coefficients  of  the  sample.  The 
frequency  domain  apparatus  used  in  our  lab  to  measure  tissue  optical  properties  is  sketched  in 
Figure  1.  Briefly,  an  oscillator  modulates  a  laser  diode  at  a  radio  frequency  (RF)  of  100  MHz. 
This  modulated  light  is  carried  to  the  sample  through  an  optical  fiber,  and  the  resulting  scattered 
light  is  collected  by  a  different  optical  fiber,  which  guides  the  signal  to  the  detector.  This  signal 
is  amplified  and  mixed  with  a  part  of  the  original  source  in  order  to  extract  the  amplitude  and  RF 
phase  of  the  modulated  light.  By  measuring  the  changes  in  the  DC  and  RF  amplitudes  and  the 
RF  phase  of  the  detected  signal  as  a  function  of  source-detector  separation,  we  can  determine  the 
absorption  and  scattering  coefficients  of  the  measured  tissue  region. 1  This  frequency-domain 
system  has  a  precision  of  better  than  1%  for  measurements  on  homogeneous  tissue  phantoms.2 
We  have  made  a  number  of  improvements  to  this  measurement  system,  which  will  be  discussed 
in  the  next  section. 
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Figure  1.  Schematic  diagram  of  the  frequency-domain  imaging  system. 


We  have  constructed  a  continuous  wave  immersion  imaging  system  to  examine  dynamic 
physiological  behavior  of  rodents.  A  schematic  of  the  instrument  used  for  immersion  imaging  is 
shown  below  in  Figure  2.  The  imaging  system  is  composed  of  a  light  source,  an  immersion  box, 
and  a  camera.  The  light  source  is  made  up  of  arrays  of  bright  light  emitting  diodes  (LED)  that 
emit  radiation  with  peak  intensities  at  three  wavelengths:  570  mn,  780  mn,  or  840  mn  (Epitex). 
This  diffuse,  CW  light  source  is  directed  at  the  sample  immersion  box,  which  contains  the 
animal  under  study  immersed3  in  a  heated  (37  -  38°  C),  matching  medium  (tissue  phantom) 
composed  of  water,  ink,  sugar,  and  submicrometer  polymer  spheres  (  Ropaque,  Rohm  and  Haas 
Company),  which  approximates  the  scattering  and  absorptive  properties  of  the  mouse  tissue. 

The  front  of  the  immersion  box  is  imaged  onto  the  camera.  Images  at  each  individual 
wavelength  are  then  collected,  digitized,  and  sent  to  the  computer  for  analysis. 
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Immersion  Box 


Figure  2.  Schematic  diagram  of  the  CW  immersion  imaging  system. 

The  780  nm  and  840  nm  wavelengths  lie  on  either  side  of  the  isosbestic  point  for 
hemoglobin  and  are  intended  to  provide  information  on  hemoglobin  content  and  hemoglobin 
oxygenation.4  The  570  nm  wavelength  is  very  strongly  attenuated  by  the  hemoglobin  in  the 
animal  tissues  and  is  much  less  attenuated  in  the  tissue  phantom.  This  wavelength  is  intended  to 
provide  a  images  showing  the  position  and  outline  of  the  animal  while  immersed. 

TASK  2.  IMPROVE  IMAGING  SYSTEMS 

We  have  made  substantial  improvements  to  both  the  frequency  domain  and  the  CW 
imaging  systems.  For  the  frequency  domain  system,  we  have  rebuilt  the  detection  electronics. 
The  detector  is  a  new,  multi-alkali,  miniaturized  photomultiplier  tube  (PMT)  from  Hamamatsu. 
This  tube  is  more  sensitive  to  the  NIR  wavelengths  used  in  our  experiments.  It  is  also  smaller  in 
size,  and  has  a  lower  electron  transit  time  spread.  The  output  from  the  PMT,  a  time-varying 
current,  is  directly  converted  and  amplified  to  a  voltage  via  a  wide-bandwidth  transimpedance 
amplifier  (Philips).  This  new  detection  scheme  eliminates  two  amplifier  stages,  filters,  and 
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attenuators  present  in  the  previous  apparatus,  and  increases  the  signal-to-noise  (S/N)  ratio  of  the 
detected  signals. 

Most  of  the  improvements  in  the  CW  system  have  been  aimed  at  complete  automation  of 
the  experimental  apparatus.  We  have  added  mass  flow  controllers  (MFC)  to  the  gas  lines,  which 
control  the  exact  flow  of  the  individual  gas  components  (air,  O2,  and  CO2)  to  the  mouse  nose 
cone.  These  MFCs  are  controlled  by  the  computer  through  a  Labview  program,  which  also 
controls  the  LEDs  and  the  image  acquisition.  With  the  current  setup,  we  can  deliver  and  monitor 
the  gas  concentrations  delivered  to  the  mouse,  monitor  the  temperature  of  the  immersion  bath, 
and  control  the  LED  and  camera  parameters  throughout  the  duration  of  the  experiment.  By 
carefully  controlling  and  monitoring  the  experimental  variables  in  our  system,  we  have  greatly 
reduced  noise  associated  with  the  physiological  response  of  the  mouse  to  its  environment.  Such 
improvements  have  been  essential  for  eliminating  external  sources  of  error  to  our  measurements. 
We  have  also  upgraded  the  NIR  light  source  to  increase  the  amount  of  light  incident  upon  the 
camera.  To  this  end,  we  have  added  3  times  the  number  of  LEDs  at  780  and  840  mn.  We  have 
also  added  a  Fresnel  lens  onto  the  front  of  the  immersion  box  to  more  effectively  couple  the 
scattered  light  transmitted  from  the  light  box  onto  the  CCD  camera.  This  has  allowed  us  to 
increase  the  acquisition  frame  rate  and  more  carefully  probe  the  physiological  dynamics  of  the 
tissue  vasculature. 

To  improve  the  sensitivity  and  speed  of  the  animal  imaging  experiments  we  have 
replaced  the  Point  Grey  Instruments  Research  Dragonfly  camera  with  a  Qlmaging  Retiga 
camera.  Between  the  increased  area  of  the  CCD  and  the  higher  near  infrared  quantum 
efficiency,  this  change  has  led  to  an  increased  sensitivity  by  a  factor  of  about  8.  The  Labview 
program  that  controls  the  experiment  was  reprogrammed  to  use  the  new  camera. 
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TASK  3.  PREPARE  AND  MAINTAIN  ANIMAL  MODELS 


The  overall  goal  of  this  research  is  to  improve  the  detection  of  breast  cancer  through  a 
non-invasive  vascular  function  imaging  system  using  animal  cancer  models.  In  vivo  xenograft 
models  are  useful  for  these  studies  because  they  allow  monitoring  development  (the  progressive 
vascular  invasion  by  tumor  cells),  including  neoplastic  contrast  on  range  stage  (based  on  tumor 
size).  These  xenograft  models  also  allow  us  to  examine  a  wide  variety  of  tumors  and  their  unique 
growing  characteristics  in  order  to  determine  the  strengths  and  weaknesses  of  the  current  non- 
invasive  imaging  technique.  As  a  model  for  breast  cancer  detection,  the  mice  with  tumor 
xenografts  are  partially  immersed  in  liquid  tissue  phantoms  (media)  that  simulate  the  optical 
properties  of  breast  tissue.  Although  this  approach  does  not  completely  include  the  properties  of 
tissue  heterogeneity  in  the  breast,  the  use  of  xenograft  models  are  the  most  practical  methods  for 
studying  contrast  over  a  range  and  verity  of  tumors.  In  the  course  of  study  we  also  included 
other  non-human  breast  cancer  cells,  to  help  improve  technique  detection  and  validation  of  the 
imaging  system  and  protocols. 

In  this  work  we  used  two  human  breast  cancer  models,  an  estrogen-dependent  cell  line 
(MCF-7)  and  estrogen-independent  cell  line  (MDA-MB-231)  and  three  other  cell  lines  (mouse 
embryonic  fibrosarcoma,  U87  brain  cancer,  and  DLD-1  colon  cancer).  Our  imaging  technique 
examines  the  vasculature  of  the  tumor  through  its  response  to  inhalation  of  carbon  dioxide  and 
oxygen.  The  additional  cell  lines,  particularly  the  mouse  embryonic  fibrosarcoma  and  brain 
cancer  cell  lines,  have  a  larger  amount  of  vasculature  than  the  two  breast  cell  lines.  In  patients, 
we  can  expect  considerable  variability  in  the  amount  of  tumor  vasculature.  Thus,  including  the 
additional  cell  lines  provides  a  better  indication  of  the  range  of  response  for  breast  cancer 
detection  in  patients. 
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We  had  originally  proposed  to  use  two  animal  models  for  these  studies:  a  human  cancer 
cell  line  in  mice  and  a  chemically  induced  tumor  in  rats.  Because  of  the  significance  of 
investigating  a  broad  range  of  tumor  vasculature,  we  have  discussed  with  the  Grants  Manager  for 
this  project  substituting  the  four  additional  cell  line  models  for  the  one  chemically-induced  tumor 
model.  The  Grants  Manager  allowed  that  this  was  a  small  enough  change  that  a  formal  change 
to  the  Statement  of  Work  was  not  required. 

TASK  4.  PERFORM  ANIMAL  IMAGING 

Imaging  experiments  were  conducted  on  animals  with  tumor  volume  of  300-1000  mm3. 
For  each  experiment  2-4  animals  were  used.  Prior  to  imaging,  the  animals  were  anesthetized 
with  40  mg/kg  pentobarbital.  The  mice  are  then  secured  to  a  3-inm  Plexiglas  platfonn  with  clear 
surgical  tape.  The  tape  used  is  an  improvement  over  the  earlier  technique,  which  employed 
black  nylon  tape  as  an  adhesive.  The  black  tape  produced  regions  of  contrast  following 
subtraction  of  a  background  image,  described  below.  Because  we  are  using  a  differential 
measurement  to  examine  changes  in  tissue  optical  properties,  any  movement  of  the  animal 
generates  artifacts  in  the  subtracted  measurements.  Therefore,  anesthesia  is  given  in  further 
doses  of  20  mg/kg  as  needed  to  reduce  stress  associated  with  immersion  and  to  keep  the  animal 
immobilized.  Using  the  new  computer-controlled  apparatus,  custom  gas  mixtures  were 
administered  to  the  immersed  mouse  via  a  nose  cone  at  a  flow  rate  of  approximately  4.2  1/min. 
Such  a  high  flow  rate  was  chosen  to  eliminate  rebreathing  artifacts.  Typically,  we  cycled 
between  pure  air,  pure  Cfi,  and  carbogen  (95%  Cfi,  5%  CCfi).  In  order  to  test  the  different  effects 
of  hyperoxia  (elevated  Cfi)  versus  those  of  hypercapnia  (elevated  CO?),  we  also  have 
experimented  with  varying  the  levels  of  CO?  in  the  presence  of  O?  and  air,  respectively.  The 
optical  path  length  of  the  immersion  box  is  adjusted  to  match  the  thickness  of  the  mouse  (~2- 
2.5  cm).  With  this  thickness,  the  exposure  time  of  the  camera  allows  us  to  measure  both  780  nm 
and  840  nm  wavelengths  at  approximately  one  frame  per  second.  At  the  end  of  the  experiment 
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the  animals  were  sacrificed  with  an  overdose  of  anesthetic  agent  followed  by  cervical 
dislocation. 

A  background  image  of  individual  mice  was  recorded  for  each  wavelength  at  the 
beginning  of  each  experiment.  Following  background  collection,  the  mouse  was  imaged  during 
and  after  the  administration  of  the  different  gas  protocols.  From  each  of  these  images  the 
background  image  was  subtracted  to  produce  a  differential  image  of  the  response  of  the  optical 
properties  of  the  mouse  to  the  different  gases.  Figure  3,  below,  demonstrates  the  enhancement  of 
the  imaging  contrast  with  this  differential  method.  The  first  panel  shows  one  of  these  static 
images  recorded  during  the  administration  of  oxygen  gas.  The  outlines  of  both  the  mouse  and 
the  tumor,  determined  by  recording  an  image  of  the  animal  without  the  immersion  medium,  have 
been  placed  on  top  of  the  image  as  a  guide  to  the  eye.  The  mouse’s  head  is  out  of  the  liquid 
immersion  medium  and  is  at  the  top  of  the  field  of  view  in  the  figure.  The  second  panel  shows 
the  same  image  after  subtraction  of  the  background.  Although  the  actual  boundaries  of  the 
mouse  and  tumor  are  obscured  by  the  scattering  of  the  immersion  medium,  it  is  clear  from  the 
difference  image  in  Figure  3  b  that  there  are  distinct  regions  of  contrast  between  the  tumor  and 
the  surrounding  tissues  of  the  mouse.  While  there  is  a  region  of  enhanced  contrast  in  the  first 
panel  due  to  the  static  optical  properties  of  the  tumor,  it  is  also  obvious  that  the  contrast  in  this 
region  is  much  greater  in  the  differential  image. 


Figure  3.  Images  of  mouse  at  840  nm  at  100  s  after  administration  of  carbogen  gas;  a)  static 
image;  b)  image  from  3a  with  background  subtracted. 
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Principal  Component  Analysis 

The  imaging  experiments  described  above  generated  large  sets  of  data.  Typically,  images 
with  I  05  pixels  at  two  wavelengths  are  recorded  every  2-10  seconds  over  the  cycling  period  of 
carbogen  administration  (approximately  10  to  20  minutes).  Based  on  these  experimental  results, 
we  expect  to  see  <  7%  change  in  image  intensity  following  carbogen  administration.  Because 
extracting  such  small  signal  changes  from  large  data  sets  poses  a  formidable  challenge, 
researchers  have  developed  techniques  that  generate  smaller  sets  of  orthogonal  images  to 
describe  the  generated  data.5'6  In  practice,  these  methods  have  been  shown  to  accurately  describe 
data  sets  of  10,000  images  with  only  -100  eigenimages. 

In  the  most  basic  adaptation  of  these  methods,  known  as  principal  component  analysis 
(PCA),  the  set  of  recorded  images  is  represented  by: 

/=/0,x)  (1) 

where  x  describes  the  spatial  pixel  grayscale  values  of  the  image,  and  1  is  the  time  at  which  the 
image  data  was  collected.  Researchers  have  shown  that  these  images,  f(t,x),  can  be  decomposed 
into  the  set  of  orthogonal  functions  a„(t)  and  cpn( x )  by: 

/(t,x)  =  2^„«„(0<P„(x).  (2) 

n 

A  series  of  T  time  images  containing  P  pixels  can  be  described  by  the  matrix: 

/(U)  /( 1,2)  ...  f(l,P) 

f  (2,1)  /  (2,2)  ...  f(2,P) 

f(T,  1)  .  f(T,P ) 

This  matrix  can  then  be  decomposed  into  the  different  a„(t)  and  <p„(x)  components  through  the 
general  technique  of  singular  value  decomposition: 
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M=AUVt  (5) 

The  columns  of  V  contain  the  orthonormal  spatial  basis  functions,  the  orthonormal  columns  of  A 
describe  the  time-dependence  of  the  spatial  basis  functions,  and  U  contains  the  weighting  factors 
for  the  two  matrixes  A  and  V. 

As  a  first  step  in  processing  the  data,  we  apply  this  simplified  PCA  method  to  determine 
changes  in  oxyhemoglobin  and  deoxyhemoglobin,  scaled  by  some  pathlength  factor  /  (see  Kotz 
et  al.  2006,  included  in  appendix).  The  time-dependent  images  that  describe  A[Hb]  and  A[HbC>2\ 
were  ordered  into  a  matrix  as  shown  in  equation  (3),  and  the  singular  value  decomposition  was 
carried  out  to  obtain  the  matrices  A,  U,  and  V.  Figure  4  presents  a  plot  of  the  normalized  scaling 
factors  contained  along  the  diagonal  of  U.  Only  the  first  three  or  four  eigenimages  contribute 
significantly  to  the  set  of  images  that  describe  the  hemoglobin  dynamics  in  our  study. 

Figure  5  shows  the  first  eigenimage  corresponding  to  the  first  two  columns  of  matrix  V. 
The  contrast  between  the  tumor  and  the  surrounding  tissue  is  evident  in  the  second  image.  The 
time-dependent  weighting  of  the  second  eigenimage  in  the  A  [ Hb](t )  and  A  [HbC>2](t)  sets  of 
images  can  be  determined  from  the  matrix  product  of  A*U,  and  is  shown  in  Figure  6. 
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Figure  4.  Normalized  scaling  factors  (eigenvalues)  for  principal  component  analysis  of  image 
stack. 


First  Eigenimage 


Figure  5.  First  eigenimage  from  principal  component  analysis. 
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Figure  6.  Time-dependent  weighting  of  the  first  eigenimage  from  principal  component  analysis. 

TASK  6.  COMPARE  DYNAMIC  AND  STATIC  MEASUREMENTS 

The  response  of  the  mouse  to  the  gas  inhalation  protocols  was  assessed  by  examining  the 
time-dependent  sequence  of  images  recorded  at  each  wavelength.  As  a  first  step  in  the  data 
analysis,  we  have  recorded  the  mean  and  standard  deviation  of  each  background-subtracted 
image  through  time.  The  variation  in  pixel  intensity  through  time  gives  a  quantitative  measure  of 
the  changes  in  the  images  in  response  to  the  inhaled  gases.  This  variation  in  pixel  intensity  is 
shown  in  Figure  7  below  for  each  of  the  wavelengths  studied.  The  response  is  clearly  enhanced 
in  the  regions  of  the  tumor,  and  it  appears  as  though  the  there  is  a  enhancement  in  the  response  of 
the  tumor  to  the  gases.  The  general  trends  resulting  from  gas  inhalation  are  consistent  with  point 
measurements  observed  by  other  researchers.7 

The  dynamics  of  the  enhanced  contrast  between  the  tumor  tissue  and  the  mouse  tissue 
due  to  the  inhalation  of  the  different  gases  was  monitored  by  spatially  averaging  the  changes  in 
the  image  intensity  over  areas  corresponding  to  signals  due  to  cancerous  tissue  and  areas  due  to 
normal  tissues  within  the  mouse.  These  averaged  data  are  shown  in  Figure  8  for  differences  in 
the  780  nrn  and  840  mn  images,  respectively.  The  dashed  line  represents  changes  in  the  tumor 
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tissue,  the  thick  solid  line  is  for  an  adjacent  region  within  the  mouse  that  does  not  contain  the 
tumor,  and  the  line  represents  the  average  of  a  part  of  the  image  not  containing  the  mouse.  It  is 
clear  from  both  figures  that  there  are  distinct  differences  between  the  dynamics  of  the  tumor 
tissue  when  compared  with  the  normal  mouse  tissue.  It  is  also  clear  that  the  response  of  the 
tissues  is  correlated  to  the  presence  of  CO2  in  the  840  mn  images,  whereas  the  modulation  of  O2 
dominates  the  response  at  780  mn. 


780  nm  Standard  Deviation  840  nm  Standard  Deviation 


Figure  7.  Standard  deviations  of  the  time-varying  pixel  absorbance  values. 


Figure  8.  Dynamic  response  of  the  averaged  signal  due  to  the  tumor  tissue  (dashed 
line)  and  the  mouse  tissue  (solid  line)  at  780  nm  and  840  nm. 
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We  have  also  measured  the  static  properties  of  the  mouse  tissue  with  the  modified 
frequency  domain  apparatus.  For  these  measurements,  an  anesthetized  mouse  was  suspended 
between  a  U-shaped  bar  of  aluminum.  The  mouse  was  held  in  place  with  a  custom  harness.  The 
properties  of  the  immersion  media  were  measured  as  previously  discussed,  and  following  these 
measurements,  the  mouse  was  scanned  between  the  source  and  detector  fibers  while  measuring 
the  amplitude  and  phase  of  the  scattered  signal.  A  sample  scan  recorded  at  780  mn  is  shown 
below  in  Figure  9.  Although  the  immersion  media  was  not  an  exact  match  to  the  properties  of 
the  mouse,  we  were  able  to  detennine  the  average  optical  properties  of  the  mouse  tissues,  and  we 
were  able  to  observe  a  region  of  the  scan  with  anomalous  scattering,  which  was  correlated  with 
the  position  of  the  tumor.  We  have  used  these  results  to  further  refine  the  exact  scattering  and 
absorption  properties  of  the  immersion  media  used  in  subsequent  measurements. 


Figure  9.  Scan  of  mouse  tissue  optical  properties  at  780  nrn  taken  with  frequency  domain  apparatus. 


17 


PERSONNEL 


The  following  personnel  were  supported  by  this  research  effort: 


Name 

Role 

Duties 

Gregory  Faris 

Principal  Investigator 

Project  Leadership 

Kenneth  Kotz 

Postdoctoral  Fellow 

Experimental  Studies 

Konstantinos  Kalogerakis 

Investigator 

Experimental  Studies 

Khalid  Amin 

Investigator 

Animal  Models 

Juan  Orduna 

Investigator 

Animal  Models 

Zhaohui  Wang 

Investigator 

Animal  Models 

Wan-Ru  Chao 

Investigator 

Animal  Models 

Barbara  Sato 

Investigator 

Animal  Models 

Harold  Javitz 

Investigator 

Statistician 

David  Watters 

Investigator 

Instrumentation  Consulting 

Steven  Young 

Technician 

Electronics  Support 

Karl  van  Dyk 

Technician 

Mechanical  Support 

William  Olson 

Technician 

Mechanical  Support 

Ross  Jackson 

Technician 

Mechanical  Support 

Gabe  Hernandez 

Administrative  Assistant 

Technical  Document  Research 

18 


STUDENTS 


The  following  students  contributed  to  this  research  effort  under  support  of  the  Research 
Experiences  for  Undergraduates  Program  of  the  National  Science  Foundation  (NSF  Physics). 

Kendra  Rand,  a  junior  at  Carthage  College,  Kenosha,  WI 

William  Boenig,  a  junior  at  Santa  Clara  University,  Santa  Clara,  CA 

Nevena  Rakuljic,  a  sophomore  at  University  of  California  San  Diego 

Ashley  Gibbs,  a  sophomore  at  Whitworth  College 


19 


KEY  RESEARCH  ACCOMPLISHMENTS 


We  have  achieved  the  following  accomplishments  on  this  project: 

•  Developed  a  simple  and  practical  cw  animal  imaging  system. 

•  Developed  computer-controlled  gas  control  system  for  easily  variable  animal 
inhalation  experiments  with  on-demand  gas  mixing. 

•  Tested  our  approach  on  5  animal  tumor  models. 

•  Compared  static  and  dynamic  contrasts  mechanisms. 

•  Developed  principal  component  analysis  tools  for  image  stack  analysis. 

•  Shown  dramatic  contrast  enhancement  for  our  dynamic  imaging  method  that  is 
the  basis  for  clinical  studies. 
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CONCLUSIONS 


We  have  shown  that  our  dynamic  imaging  method  can  provide  strong  cancer-specific 
contrast  in  animal  tumor  models.  Based  on  this  project,  funding  is  now  in  place  for  clinical 
studies  of  this  non-invasive  imaging  method.  If  functional  optical  imaging  can  enhance  early 
detection  of  breast  cancer  by  detecting  tumors  missed  on  x-ray  mammography  (false  negative 
results),  it  could  reduce  the  economic  and  human  cost  associated  with  later  detection  of  disease. 
If  this  approach  can  improve  the  diagnosis  of  breast  cancer  (reduce  the  number  of  false  positive 
diagnoses),  it  could  reduce  the  worry  and  economic  cost  of  unnecessary  biopsies.  By  providing 
more  specific  functional  image  information,  optical  imaging  can  enhance  the  ability  to  detect  and 
diagnose  cancer,  compared  with  what  is  achievable  with  the  physical  image  information 
provided  by  x-rays.  For  example,  because  of  the  low  cost  of  optical  instrumentation,  optical 
imaging  might  be  used  in  combination  with  x-ray  mammography,  which  should  provide  greater 
sensitivity  and  specificity  than  x-rays  alone.  With  the  transition  to  digital  x-ray  mammography, 
optical  imaging  and  x-ray  mammography  can  even  share  the  same  camera,  thereby  providing 
excellent  registration  of  the  two  modalities. 
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ABSTRACT 

We  are  using  rodent  animal  models  to  study  and  compare  contrast  mechanisms  for  detection  of  breast  cancer.  These 
measurements  are  performed  with  the  animals  immersed  in  a  matching  scattering  medium.  The  matching  scattering 
medium  or  liquid  tissue  phantom  comprises  a  mixture  of  Ropaque  (hollow  acrylic/styrene  microspheres)  and  ink.  We 
have  previously  applied  matched  imaging  to  imaging  in  humans.  Surrounding  the  imaged  region  with  a  matched  tissue 
phantom  compensates  for  variations  in  tissue  thickness  and  geometry,  provides  more  uniform  illumination,  and  allows 
better  use  of  the  dynamic  range  of  the  imaging  system.  If  the  match  is  good,  the  boundaries  of  the  imaged  region  should 
almost  vanish,  enhancing  the  contrast  from  internal  structure  as  compared  to  contrast  from  the  boundaries  and  surface 
topography.  For  our  measurements  in  animals,  the  immersion  plays  two  additional  roles.  First,  we  can  readily  study 
tumors  through  tissue  thickness  similar  to  that  of  a  human  breast.  Although  the  heterogeneity  of  the  breast  is  lost,  this  is 
a  practical  method  to  study  the  detection  of  small  tumors  and  monitor  changes  as  they  grow.  Second,  the  immersion 
enhances  our  ability  to  quantify  the  contrast  mechanisms  for  peripheral  tumors  on  the  animal  because  the  boundary 
effects  on  photon  migration  are  eliminated.  We  are  currently  developing  two  systems  for  these  measurements.  One  is  a 
continuous-wave  (CW)  system  based  on  near-infrared  LED  illumination  and  a  CCD  (charge-coupled  device)  camera. 
The  second  system,  a  frequency  domain  system,  can  help  quantify  the  changes  observed  with  the  CW  system. 

Keywords:  Compensated  transillumination,  frequency-domain  measurements,  animal  models,  breast  cancer,  tumor 
detection,  in  vivo  imaging. 


1.  INTRODUCTION 

The  observation  of  illumination  through  tissue,  known  as  transillumination  or  diaphanography,  was  one  of  the  first 
methods  employed  for  medical  imaging,  dating  back  to  the  early  19th  century.1’2  Optical  mammography  was  first 
performed  in  1929. 3  The  introduction  of  cooler  light  sources  and  the  use  of  film  or  a  camera  to  detect  near-infrared 
wavelengths  generated  considerable  research  activity  in  breast  transillumination  during  the  1970’s  and  1980’s.4’6 
Nevertheless,  optical  mammography  proved  inferior  to  x-ray  mammography.7-8  The  primary  problem  with  optical 
mammography  is  poor  spatial  resolution,  which  reduces  contrast  in  smaller  tumors.  However,  this  limitation  can  be 
reduced  or  even  overcome  by  providing  functional  (physiological)  imaging  information. 

There  has  been  a  resurgence  of  interest  in  optical  techniques  for  breast  cancer  detection  in  the  past  decade,  spurred  on  by 
better  understanding  of  light  propagation  through  tissue,9-10  time-resolved  optical  techniques  for  imaging,11  15  better 
quantification  of  tissue  properties,16-17  and  new  image  reconstruction  algorithms.18'20  This  interest  has  been  reflected  in 
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the  development  of  industrial  prototype  imaging  systems  at  Philips, 12  Zeiss, 14  and  Siemens. 15  Applications  of  advanced 
optical  mammography  techniques  have  shown  promising  results.14’21'24  Several  research  groups  are  examining  the  use 
of  contrast  agents  for  breast  cancer  detection.24'28  Others  have  performed  animal  studies  using  window  chambers  to 
provide  direct  optical  access  to  study  the  microenvironment  of  a  tumor.29’30  A  number  of  ex  vivo  studies  of  the  optical 
properties  of  normal  and  diseased  breast  tissue  have  also  been  performed.31'33 

One  limitation  for  in  vivo  optical  imaging  is  due  to  non-uniform  thickness  and  variable  geometry.  Small  variations  in 
the  thickness  of  the  illuminated  tissue  lead  to  large  variations  in  the  intensity  of  transmitted  light.  Asymmetries  in  the 
location  and  orientation  of  the  tissue  surface  cause  variations  in  the  reflectivity  of  the  tissue.  Additionally,  light  may  be 
lost  or  diffuse  away  at  tissue  boundaries  transverse  to  the  imaging  plane.  In  combination,  these  effects  cause  significant 
variations  in  the  intensity  of  the  transmitted  light  and  reduce  contrast  within  the  observed  tissue.  In  some  cases,  the 
undesired  effects  of  intensity  variation  can  be  reduced  by  measuring  the  optical  path  length  and  by  making  the 
corresponding  adjustments  to  the  image  contrast.  This  approach  has  been  applied  in  frequency-domain  techniques 
where  the  phase  shift  of  an  amplitude-modulated  light  beam  was  used  to  determine  the  optical  path  length.34  In  our 
experiments,  the  uniformity  and  contrast  of  the  image  is  improved  by  immersing  the  imaged  tissue  in  a  liquid  medium 
with  the  same  absorption  and  scattering  behavior  as  the  tissue  under  study.  Such  matching  of  the  optical  properties  by 
immersion  results  in  a  more  uniform  image  and  in  optimized  contrast.  When  a  good  match  is  achieved,  the  tissue 
boundaries  will  practically  disappear,  and  the  contrast  in  the  image  will  mainly  reflect  variations  in  internal  structure.  In 
addition,  imaging  in  a  matching  liquid  medium  usually  results  in  improved  and  more  uniform  illumination  conditions. 

Matching  techniques  have  been  used  in  the  past  for  the  determination  of  tissue  properties16’35’36  and  for  scanned  time- 
resolved  optical  mammography.37’38  Previous  work  at  SRI  International  has  combined  CW  imaging  with  immersion  in  a 
matching  medium  to  obtain  images  of  the  human  hand.39  A  conventional  incandescent  lamp  illuminated  a  Plexiglas  box 
containing  the  matching  medium,  a  mixture  of  Ropaque,  water  and  methylene  blue  dye.  Matching  between  the  human 
tissue  and  relative  concentration  of  the  dye  and  Ropaque  was  achieved  with  frequency-domain  measurements.36’40 
Images  of  a  human  hand  immersed  in  the  matching  medium  were  recorded  with  an  intensified  CCD  and  a  conventional 
camera.  Comparisons  with  images  recorded  without  immersion  showed  that  images  with  immersion  display  more  detail 
and  internal  structure  than  can  be  seen  without  the  matching  conditions.39 

While  x-ray  imaging  primarily  reveals  structural  information,  optical  imaging  provides  both  structural  information  and 
information  on  tissue  function  through  the  use  of  optical  spectroscopy.  Spectroscopic  information  may  be  enhanced 
through  application  of  time-resolved  methods  such  as  time-domain  or  frequency-domain  techniques  to  provide 
quantitative  information  on  chemical  state  or  composition.  For  example,  optical  measurements  performed  at  different 
wavelengths  can  provide  information  on  total  hemoglobin  content  and  hemoglobin  oxygenation.42  This  type  of 
functional  information  is  significant  for  breast  cancer  detection.  The  angiogenesis  associated  with  tumor  growth 
typically  leads  to  elevated  local  hemoglobin  concentrations.22,23’43  In  addition,  tumors  exhibit  atypical  vasculature  and 
vasoactivity,  and  are  often  hypoxic,  a  condition  that  can  be  observed  optically  by  hemoglobin  oxygenation 
measurements.  Tumor  hypoxia  may  also  be  used  to  guide  treatment  because  more  hypoxic  tumors  are  likely  to  be 
metastatic  or  invasive.44'46  The  distinct  morphology  of  tumors  has  also  been  shown  to  provide  a  source  of  contrast 
through  variations  in  the  optical  scattering  coefficient.47 

In  this  proceeding,  we  present  results  from  preliminary  experiments  inspired  by  a  new  form  of  functional  optical 
imaging.  The  main  objective  is  to  study  the  contrast  enhancement  available  using  inspiration  of  varying  amounts  of 
carbon  dioxide  and  oxygen.  These  studies  are  being  performed  in  rodents.  The  animal  model  employed  allows  the 
study  of  contrast  during  tumor  progression  and  over  a  range  of  tumor  sizes.  We  use  imaging  by  immersion  in  a 
matching  medium  to  improve  quantification,  simplify  measurements  on  the  small  animals,  and  simulate  tissue  depth 
similar  to  that  of  the  human  breast.  These  studies  attempt  to  use  the  full  range  of  available  optical  contrast  by 
augmenting  functional  optical  imaging  using  dynamic  measurements  related  to  tumor  vascular  function  associated  with 
inspired  gases.  The  broadest  use  of  available  contrast  will  be  most  effective  for  improving  sensitivity  and  specificity  in 
tumor  detection  and  monitoring. 
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2.  EXPERIMENTAL  APPARATUS  AND  RESULTS 


We  are  using  two  optical  systems  for  our  studies — a  frequency  domain  system  for  quantification  and  a  CW  system  for 
monitoring  rapid  changes.  Both  systems  use  immersion  of  the  animals  in  liquid  tissue  phantoms  to  improve  the 
uniformity  and  contrast  of  the  images  and  facilitate  quantification  for  the  small  volumes  imaged. 

The  first  instrument  is  based  on  the  frequency  domain  approach,  a  rather  simple  and  accurate  method  for  quantitative 
optical  measurements  of  tissue  properties.  This  technique  uses  measurements  of  the  optical  delay  of  light  through  a 
medium  to  quantitatively  determine  absorption  and  scattering  properties  on  human  or  animal  tissue.  The  phase  shift  of 
the  modulated  light  increases  with  the  amount  of  scatter,  while  its  amplitude  decreases  with  absorption  by  the  medium. 
A  simplified  schematic  diagram  of  a  frequency  domain  measurement  system  developed  in  our  laboratory  is  shown  in 
Figure  1.  An  oscillator  modulates  a  light  source  at  a  radio  frequency  (RF)  of  100  MFIz.  This  light  is  carried  to  the 
measured  region  through  an  optical  fiber.  A  second  fiber  carries  scattered  light  to  a  detector.  A  mixer  detects  the 
amplitude  and  RF  phase  of  the  modulated  light.  By  measuring  the  changes  in  the  DC  and  RF  amplitudes  and  the  RF 
phase  of  the  detected  signal  as  a  function  of  source-detector  separation,  we  can  determine  the  absorption  and  scattering 
coefficients  of  the  measured  tissue  region.48  This  frequency-domain  system  has  a  precision  of  better  than  1%  for 
measurements  on  homogeneous  tiss  ue  phantoms.40 


Figure  1 .  Schematic  diagram  of  the  frequency-domain  imaging  system. 
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Once  the  absorption  and  scattering  behavior  of  the  immersion  medium  has  been  characterized,  the  only  measurements 
required  to  measure  a  tissue  sample  are  the  differences  in  phase  and  amplitude  with  and  without  the  sample  between  the 
optical  fiber  bundles.  Thus,  the  determination  of  the  absorption  coefficient,  pa,  and  the  reduced  scattering  coefficient, 
/js',  as  a  function  of  wavelength,  time,  or  position,  is  a  straightforward  matter. 

We  use  the  CW  instrument  to  examine  the  dynamic  physiological  behavior.  These  measurements  are  performed  using  a 
CW  source  and  a  two-dimensional  detector,  such  as  a  CCD  or  conventional  camera.  The  instrumentation  for  this 
immersion-based  imaging,  or  compensated  transillumination,  is  shown  in  Figure  2.  The  system  comprises  an  LED 
(light-emitting  diode)  array  as  the  illumination  source,  filters,  a  box  for  the  immersion  medium,  and  a  CCD  camera. 
Two  types  of  LED  arrays  are  currently  used  with  output  radiation  centered  at  780  and  840  nm,  respectively.  The 
imaging  system  is  suitable  for  a  fast  frame  acquisition  rate  and  averaging.  Switching  between  LED  arrays  enables 
measurements  at  different  wavelengths  and  the  determination  of  hemoglobin  content  and  hemoglobin  oxygenation.42 
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Figure  2.  Schematic  diagram  of  the  CW  immersion  imaging  system. 


The  studies  of  vascular  function  imaging  are  performed  on  animal  breast  cancer  models  (rats  with  chemically  induced 
tumors  and  nude  mice  with  tumor  xenografts).  We  use  animal  models  because  they  allow  monitoring  contrast  on  a 
range  of  tumor  sizes  and  stages  of  development.  Furthermore,  the  measurements  are  noninvasive  and  thus  can  be  readily 
repeated  on  animals  as  our  instrumentation  and  methods  are  refined.  Measurement  through  tissue  thickness  similar  to 
the  human  breast  can  be  performed  by  partially  immersing  the  anesthetized  animals  in  liquid  tissue  phantoms  that 
simulate  the  optical  properties  of  breast  tissue.  Although  this  approach  does  not  include  the  effects  of  tissue 
heterogeneity  in  the  breast,  it  is  the  most  practical  method  for  studying  contrast  over  a  range  of  tumor  sizes. 

In  this  paper,  we  present  preliminary  results  from  our  imaging  studies.  Figure  3  shows  an  image  of  a  nude  mouse 
recorded  with  our  CW  imaging  apparatus  operating  at  780  nm.  The  anesthetized  mouse  is  attached  by  its  legs  onto  a 
3mm  Plexiglass  platform  with  black  tape.  The  Plexiglass  platform  is  placed  in  the  immersionbox  in  a  manner  such  that 
the  mouse  is  illuminated  from  the  front  and  its  abdomen  is  facing  the  light  source.  The  immersion  box  is  filled  with  a 
tissue  phantom  approximating  the  absorption  and  scattering  properties  of  the  mouse.  The  tumor,  which  has  a  diameter 
of  approximately  1  cm,  can  be  distinguished  clearly  as  the  bright  circle  on  the  lower  part  of  the  left  side  of  the  image. 
The  head  of  the  mouse  is  outside  of  the  immersion  box,  above  the  surface  of  the  tissue  phantom  and  is  not  included  in 
the  image.  Part  of  the  tail  of  the  mouse  can  be  seen  at  the  bottom  of  the  image.  The  mouse  has  been  anesthetized  by 
injection  of  sodium  pentobarbital  and  is  breathing  ambient  air.  The  image  presented  is  a  single  still  frame  and  no 
averaging  has  been  performed. 
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Figure  3.  Image  of  a  nude  mouse  with  tumor  obtained  by  compensated  transillumination  at  780  nm. 


Figure  4  presents  an  intensity  profile  obtained  by  integrating  the  light  intensity  in  the  rectangular  area  shown  in  Fig.  3. 
The  central  part  of  the  mouse’s  body  is  the  darker  region  and  the  tumor  the  brighter  one.  A  slope  exists  in  the  baseline 
illumination  intensity  between  the  left  and  right  part  of  the  image.  Future  experiments  will  address  this  issue  and  include 
improvements  in  the  instrumentation  so  as  to  achieve  better  illumination  uniformity  and  image  symmetry.  Nevertheless, 
even  from  a  single  frame,  it  is  possible  to  identify  readily  the  location  and  size  of  the  tumor  as  well  as  the  relative 
hemoglobin  concentration. 


Figure  4.  Caption  for  profile  graph  with  cross  section  through  tumor 


Proc.  of  SPIE  Vol.  4955 


149 


The  appearance  of  the  tumor  in  Fig.  3  as  a  bright  spot  indicates  lack  of  blood  content  relative  to  the  surrounding  tissue. 
This  result  is  most  likely  a  manifestation  of  tumor  necrosis  and  indicates  that  the  imaging  is  capable  of  distinguishing 
between  different  types  of  tumors,  e.g.,  a  tumor  deprived  of  blood  vs.  one  filled  with  blood.  Our  apparatus  allows  the 
study  of  the  progression  of  the  tumor  and  monitoring  of  any  necrotic  tissue  as  it  expands.  In  future  experiments,  we  plan 
to  investigate  these  possibilities  by  studying  different  types  of  tumors,  by  performing  an  autopsy  on  the  animal  after  the 
experiments,  by  using  more  than  one  wavelength  for  the  illumination,  and  by  administering  carbogen  gas  to  the  animal 
under  study. 


4.  DISCUSSION 

The  vasculature  produced  through  angiogenesis  in  tumors  has  atypical  characteristics,  which  in  turn  provide  the  basis  for 
our  imaging  approach.  Blood  vessels  in  tumors  are  often  highly  abnormal,  exhibiting  distended  capillaries  with  leaky 
walls  and  sluggish  flow.  These  properties  of  tumor  vasculature  provide  at  least  three  different  types  of  contrast  for 
optical  imaging  in  conjunction  with  varying  levels  of  inspired  oxygen  and  carbon  dioxide.  These  types  of  contrast  are 
due  to  low  oxygenation  improvement,  atypical  vasoactivity,  and  blood  pooling.  First,  when  elevated  levels  of  oxygen  are 
administered,  the  oxygenation  rise  in  tumors  is  typically  much  lower  than  that  of  normal  tissue  types.50'52  Second,  tumor 
vasculature  often  shows  vasoactivity  in  opposition  to  that  of  normal  tissues,  for  example  exhibiting  vasodilation  while 
normal  tissue  vasoconstricts.53’54  It  has  been  shown  that  tumor  vasculature  vasoconstricts  in  response  to  carbogen  (95% 
02/5%  COi)  administration,55’56  which  is  the  opposite  of  what  would  be  anticipated.  Further  evidence  of  compromised 
microvascular  control  has  also  been  observed  through  measurements  of  NADH  signals  by  Dr.  Britton  Chance  and 
collaborators.57  As  both  oxygen  and  carbon  dioxide  are  vasoactive,  atypical  tumor  behavior  from  administration  of 
changing  levels  of  these  gases  should  provide  strong  imaging  contrast.  Finally,  as  tumor  vessels  are  often  contorted  and 
leaky,  any  blood  pooling  in  these  vessels  will  lead  to  delayed  response  to  oxygenation  changes,  providing  another  good 
contrast  mechanism. 

Additional  facts  motivate  investigations  of  dynamic  contrast  associated  with  tumor  vascular  function.  First,  the 
breast  is  highly  heterogeneous,  comprising  the  lobes  (glandular  tissue),  fat,  connective  tissue,  ducts,  and  the  supporting 
vasculature.  We  believe  that  using  a  broader  palette  of  contrast  mechanisms  will  provide  more  specificity  for  optical 
imaging  and  help  compensate  for  the  heterogeneity  of  the  breast.  Second,  the  more  successful  noninvasive  optical 
measurements,  such  as  pulse  oximetry  and  functional  brain  imaging,  are  dynamic  in  nature.  In  pulse  oximetry,58’59  the 
dynamic  signal  of  the  cardiac  pulse  is  combined  with  the  optical  signals  of  oxygenated  and  deoxygenated  hemoglobin  to 
derive  blood  oxygenation  levels  without  detailed  knowledge  of  the  tissue  type  and  level  of  scattering.  For  brain 
imaging,60'63  dynamic  measurements  of  hemoglobin  or  hemoglobin  oxygenation  are  performed  during  specific  mental 
tasks  such  as  finger  tapping  or  visual  stimulus.  The  successes  of  these  approaches,  which  combine  dynamic  and 
functional  information,  lead  us  to  apply  a  similar  combination  of  dynamic  and  functional  measurements  for  the  studies 
proposed  here.  Finally,  recent  theoretical  work  has  demonstrated  improved  results  from  differential  optical  imaging64’65 
Dynamic  and  differential  imaging  techniques  are  similar  in  that  they  both  rely  on  changes  in  optical  contrast  over  time. 

By  providing  more  specific  functional  image  information,  we  believe  that  optical  imaging  can  enhance  the  ability  to 
detect  and  diagnose  cancer,  compared  to  what  is  achievable  with  the  physical  image  information  provided  by  x-rays.  In 
addition,  because  of  the  low  cost  of  the  optical  instrumentation  we  are  using,  optical  imaging  could  be  used  easily  in 
combination  with  x-ray  mammography,  and  their  synergy  should  provide  greater  sensitivity  and  specificity  than  x-rays 
alone. 
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ABSTRACT 


We  have  studied  the  dynamic  changes  in  tissue  vasculature  following  the  inhalation  of  hyperoxic  gasses  in 
rodents  as  a  model  for  optical  breast  cancer  detection.  We  have  used  a  CW  apparatus  to  measure  the  near 
infrared  (NIR)  optical  properties  of  animal  models  immersed  in  a  liquid  tissue  phantom.  By  looking  at  the 
transmission  of  different  wavelengths  in  the  NIR,  we  were  able  to  qualitatively  observe  changes  in  blood 
oxygenation  following  the  inhalation  of  different  mixtures  of  C02  and  02.  These  changes  enhanced  the  image 
contrast  between  cancerous  tissue  and  normal  tissues  of  the  rodents.  The  oxygenation  dynamics  of  the  tumors, 
following  inhalation  of  the  hyperoxic  gases,  exhibited  differences  from  surrounding  tissues  in  both  the 
magnitude  of  the  observed  signal  change  and  the  dynamic  response. 

Keywords:  Compensated  transillumination,  animal  models,  tumor  vasculature,  breast  cancer,  tumor  detection, 
in  vivo  imaging  contrast  enhancement. 


1.  INTRODUCTION 

The  near  infrared  (NIR)  region  between  700  nm  and  900  nm  provides  a  relatively  transparent  spectral  window 
for  observing  physiologically  important  chemical  compounds.  Two  of  the  important  NIR  absorbers  in  tissue  are 
oxygenated  ( Hb02 )  and  deoxygenated  hemoglobin  (Hb)  and  thus  this  region  of  the  spectrum  is  ideal  for 
measuring  tissue  oxygenation.1  It  is  well  established  that  tumors  exhibit  abnormalities  in  vasculature  such  as 
abrupt  endings  of  vessels,  abnormal  branching,  and  increased  leakiness  (permeability)."  Many  tumors  are  also 
often  hypoxic  which  is  a  powerful  stimulus  for  the  angiogenesis  (development  of  new  blood  vessels)  leading  to 
growth  and  invasiveness  of  the  tumor.  Many  researchers  have  attempted  to  improve  the  oxygenation  of 
cancerous  tissues  through  the  inhalation  of  increased  levels  of  oxygen,  which  makes  them  more  vulnerable  to 
certain  treatments  such  as  radiation.3  Tissue  oxygenation  can  be  measured  with  a  number  of  standard 
techniques,  however,  results  are  typically  referenced  to  an  oxygen  electrode  invasively  embedded  within  the 
tumor.4  Recently,  Liu  et  al.  have  attempted  to  monitor  such  efforts  toward  oxygenation  improvement  through 
the  in  vivo  NIR  absorption  and  p02  electrode  measurements  of  cancerous  tumors  in  rat  models.5 

Tumors  are  known  to  have  distinct  vasoactivity  in  the  presence  of  vasoactive  compounds.  Recent  studies 
have  shown  that  increased  levels  of  inspired  02  leads  to  vasoconstriction  of  all  blood  vessels,  while  increased 
levels  of  inspired  C02  leads  to  constriction  of  mature  arteries  only,  leaving  the  intratumoral  neovasculature 
unaltered.6  We  believe  that  these  properties  may  lead  to  a  mechanism  for  optical  contrast  in  static  diffuse  NIR 
transmission  images.  In  vivo  optical  imaging  of  hemodynamics  in  the  NIR  has  the  potential  to  give  not  only 
structural  information,  but  also  information  regarding  tissue  function,  as  inferred  through  the  quantitative 
information  on  the  chemical  composition  of  tissues  determined  by  optical  spectroscopy.  Thus,  we  have  recorded 
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dynamic  measurements  of  the  blood  oxygenation  parameters  of  mice  in  response  to  varying  levels  of  inspiratory 
hypercapnia  (elevated  CCb)  and  hyperoxia  (elevated  O2)  not  only  as  a  guide  to  treatment  but  also  as  a  novel 
imaging  modality  for  cancer  detection. 


2.  EXPERIMENTAL 


2.1  Animal  model 

In  this  study,  female  athymic  nude  mice  (20-25g,  Harlan  Laboratories)  were  injected  with  mouse  embryonic 
fibrosarcoma  (MEF),  U87  (Brain  Tumor)  or  human  breast  cancer  (MDA  231 )  cells  administered  subcutaneously 
on  the  dorsum  of  the  mice.  All  of  these  cell  lines  were  used  at  a  concentration  of  2-3  million  cells  in  100  pi  of 
medium  for  each  animal.  The  tumor  volumes  were  measured  twice  weekly  once  they  attain  a  measurable  size. 
In  principle,  these  xenograft  models  allow  the  study  of  contrast  during  tumor  progression  and  over  a  range  of 
tumor  sizes.  Animals  were  anesthetized  with  pentobarbital  (30-40mg/kg)  injected  intraperitoneally  at  the 
beginning  of  the  experiment. 

2.2  Imaging  apparatus 

The  instrument  used  to  measure  the  optical  properties  of  these  mice  is  shown  below  in  Fig.  1.  The 
anesthetized  animals  were  partially  immersed  in  a  temperature  controlled  (37-38°  C)  medium  that  matched  the 
optical  properties  of  the  mice  (Ropaque  and  water).7  Infrared  light  from  the  LEDs  (780  nm  and  840  nm)  was 
passed  through  the  immersion  box  and  collected  with  a  digital  monochrome  CCD  camera  (Point  Grey  Research) 
through  a  16  mm  f/0.95  lens  (JML  Optics).  Mass  flow  controllers  (MDC)  regulate  the  levels  of  oxygen  and 
carbon  dioxide  with  respect  to  ambient  air  delivered  to  the  mice  through  a  custom  nose  cone  at  constant  rate  of 

4.2  L/min.  Other  experiments  have  shown  no  effects  from  rebreathing  at  this  flow  rate.6  The  camera,  LEDs,  gas 
concentrations,  and  temperature  measurement  are  interfaced  to  and  controlled  by  a  laptop  computer.  Image  data 
is  stored  and  post-processed  using  MATLAB  and  ImageJ. 


Immersion  Box 


Fig.  1.  Experimental  Apparatus 
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3.  RESULTS 


We  used  this  apparatus  to  image  the  CW  NIR  optical  transmission  properties  of  the  animals.  The 
transmission  through  the  mice  and  the  immersion  medium  were  measured  at  the  different  wavelengths 
mentioned  above  (780  nm  and  840  nm),  while  the  gas  concentrations  were  changed.  We  cycled  between 
ambient  air,  carbogen  (95%  CF  and  5%  CCF),  100%  CF,  and  95%  air/5%  CCF-  Differential  absorption  images 
were  calculated  by  a  procedure  similar  to  that  found  in  the  literature.1'3  Briefly,  we  subtracting  the  logarithm  of  a 
baseline  (pre-protocol)  image  from  the  logarithm  the  images  at  each  time  point.  Then,  by  using  the  extinction 
coefficients  of  Hb  and  Hb02,  we  converted  the  differential  wavelength  images  into  qualitative  changes  in  the 
concentrations  of  oxygenated  and  deoxygenated  hemoglobin  (AfHb]  and  AfHbCF],  respectively).  For  these 
experiments,  we  did  not  correct  for  the  scattering  in  the  absorption  measurements,  and  thus  the  numeric  values 
given  in  the  figures  are  qualitative.  In  order  to  quickly  assess  the  overall  changes  that  occur  during  the 
experimental  run,  we  take  the  average  and  standard  deviation  of  the  pixels  through  time  corresponding  to  the 
AfHb]  and  AfFIbCF]  images.  These  are  shown  below  in  Fig  2. 


Fig.2.  Response  of  calculated  A[Hb]  and  A[Hb02]  to  gas  inhalation.  The  left  hand  images  correspond  to  A[Hb]  and 
the  right  to  AfHbCF].  The  top  images  are  the  time  averaged  mean  values  and  the  bottom  images  represent  the  standard 
deviations  from  the  mean.  The  approximate  outline  of  the  mouse  and  the  tumor  is  overlaid  on  the  images  as  a  guide  to 
the  reader.  These  were  determined  after  the  experiment  by  draining  the  tank  and  taking  an  image  that  was  registered 
to  the  dynamic  measurements. 

The  images  above  show  distinct  regions  of  high  contrast  in  the  tumor  region.  By  looking  through  time  at  the 
average  values  of  the  pixels  above  a  threshold  in  the  standard  deviation  images,  we  can  get  an  idea  of  the 
dynamics  of  the  tissue  vasculature  for  either  the  mouse  body  or  the  region  of  the  tumor  (Fig.  3).  It  can  be  seen  in 
Fig.  3  that  the  tumor  region  exhibits  a  larger  response  to  the  mouse  tissue  and  that  both  responses  are  highly 
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correlated  with  the  delivery  of  oxygen  to  the  animal  (r~0.8).  The  signal  in  both  cases  is  characterized  by  a  very 
fast  response  (<10  s)  to  the  introduction  of  the  CT  followed  by  a  very  slow  decay  back  to  the  baseline  (-30-45  s) 
These  results  are  consistent  with  similar  single-point  measurements  found  in  the  literature.3'8  Furthermore,  the 
dynamics  of  the  tumor  region  do  not  decay  back  to  the  baseline,  but  seem  to  dip  below  it.  These  are  all  general 
trends  seen  in  all  the  mice  that  we  have  studied  thus  far. 


0 


0  200 


Fig. 3.  AHb02  dynamics  in  two  different  regions  of  the  images.  The  dashed  line  represents  the  average  of  the  pixels 
that  span  the  region  of  the  lower  right-hand  image  of  Fig.  1  above  a  thr  eshold  of  0.04,  corresponding  to  changes  in  the 
tumor  tissue.  The  solid  line  represents  the  average  of  the  pixels  that  span  the  region  of  the  lower  right-hand  image  of 
Fig.  1  between  the  values  of  0.009  and  0.04,  corresponding  to  changes  in  the  mouse  tissue.  The  bar  with  light  gray 
cross-hatches  corresponds  to  100%  O2  delivery  to  the  mouse,  the  dark  gray  bar  corresponds  to  95%  O2  +  5%  CO2  gas 
delivery  to  the  mouse. 


Fig.4.  Response  of  mouse  vasculature  to  gas  inhalation  as  seen  in  the  standard  deviation  of  Hb02  difference  images. 
The  oval  in  the  image  represents  the  approximate  location  of  the  cancerous  tumor. 
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Fig.  4  shows  another  image  of  the  standard  deviations  of  pixels  for  the  AfFlbCF]  values  for  a  different 
mouse  with  a  MEF  tumor.  Again,  the  inhalation  of  CF  leads  to  the  enhancement  of  the  image  contrast  in  the 
region  of  the  tumor.  We  are  currently  working  on  the  quantifying  the  changes  in  the  AfFIbCF]  and  A|Hb] 
images,  and  on  extracting  the  dynamics  from  tumor  irrespective  of  the  surrounding  mouse  tissue.  The  current 
results,  however,  suggest  that  the  inhalation  of  gases  can  act  as  a  significant  contrast  mechanism  in  functional 
optical  imaging.  This  project  was  funded  by  the  DOD  Breast  Cancer  Research  Program,  DAMD 170920 1095 70. 
We  acknowledge  helpful  conversations  with  Dr.  Mark  Dewhirst  and  Dr.  Zishan  Haroon  of  Duke  University 
Medical  Center. 
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Inspiratory  Contrast  for  In  Vivo  Optical  Imaging 
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We  show  that  inhalation  of  oxygen  and  carbon  dioxide  leads  to  significant  contrast  for  in 
vivo  optical  imaging.  Differential  imaging  in  a  rodent  cancer  model  provides  up  to  a  factor 
of  two  variation  in  signal  change  between  tumor  and  normal  tissue.  This  exogenous 
enhancement  of  the  endogenous  contrast  due  to  oxy-  and  deoxyhemoglobin  appears  quite 
promising  for  cancer  detection. 
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OCIS  codes:  (170.0110)  Imaging  Systems;  (170.6510)  Spectroscopy,  tissue  diagnostics; 
(170.5280)  Photon  migration 

Differential  or  dynamic  optical  imaging  ’  provides  a  number  of  advantages  including  a  reduced 
sensitivity  to  boundary  effects  and  improved  sensitivity.  These  advantages  are  borne  out  in  two 
examples  of  successful  optical  measurements  in  vivo.  Functional  optical  imaging  of  the  brain3’4 
is  only  feasible  through  the  use  of  differential  imaging.  Pulse  oximetry,5’6  perhaps  the  most 
successful  clinical  use  of  light  in  human  tissue,  is  also  based  on  differential  measurements. 

During  the  process  of  angiogenesis,  tumors  develop  abnormal  vasculature,  and  as  a 
result,  cancerous  tissue  is  often  hypoxic,  a  condition  that  can  be  observed  with  hemoglobin 
oxygenation  measurements.  Hypoxic  tumors  are  also  likely  to  be  metastatic  or  invasive,  and  as  a 
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consequence,  measurements  to  monitor  blood  oxygenation  parameters  have  been  carried  out  in 
order  to  guide  the  treatment  of  such  tumors.  '  We  believe  that  it  should  be  possible  to  use  the 
endogenous  contrast  afforded  by  spectroscopic  properties  of  hemoglobin  with  the  addition  of 
exogenous  vasoactive  agents  as  a  method  to  improve  detection  of  cancerous  tumors  with  optical 
techniques. 

For  this  study,  the  differential  imaging  of  vasoactivity  is  performed  through  the 
transillumination  measurements  on  a  tumor  model  in  nude  mice.  For  the  tumor  model,  female 
athymic  nude  mice  (20-25g,  Harlan  Laboratories)  were  injected  with  mouse  embryonic 
fibrosarcoma  (MEF),  U87  (brain  tumor),  A549  (lung  tumor)  or  human  breast  cancer  (MDA  231) 
cells  administered  subcutaneously  on  the  dorsum  of  the  mice.  All  of  these  cell  lines  were  used  at 
a  concentration  of  2-3  million  cells  in  100  pi  of  medium  for  each  animal.  The  tumor  volumes 
were  measured  twice  weekly  once  they  attain  a  measurable  size.  Animals  were  anesthetized 
with  pentobarbital  (30-40  mg/kg)  injected  intraperitoneally  at  the  beginning  of  the  experiment. 
A  schematic  of  the  instrument  used  to  measure  the  transillumination  of  these  mice  is  shown  in 
Fig.  1.  The  anesthetized  animals  were  partially  immersed  in  a  temperature  controlled  (37-38°  C) 
medium  that  matched  the  optical  properties  of  the  mice  (Ropaque  and  water,  pa  =  -08  cm1,  ps  = 
12  cm’1).10  Infrared  light  from  the  LEDs  (780  nm  and  840  nm,  Epitex,  Inc.)  passed  through  the 
immersion  box  and  collected  with  a  digital  monochrome  CCD  camera  (Dragonfly,  Point  Grey 
Research)  through  a  16  mm  f/0.95  lens  (JML  Optics).  Mass  flow  controllers  (MDC)  regulate  the 
levels  of  air,  oxygen,  and  carbon  dioxide,  which  are  delivered  to  the  mice  via  a  custom  nose  cone 
at  a  constant  total  flow  rate  of  4.2  L/min.  The  camera,  LEDs,  gas  concentrations,  and 
temperature  measurement  are  interfaced  to  and  controlled  by  a  laptop  computer  on  a  mobile 
platform.  Image  data  is  stored  and  post -processed  using  MATLAB  and  ImageJ. 
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Using  this  apparatus,  CW  transillumination  measurements  were  recorded  at  each 
wavelength.  Differential  absorption  images  were  obtained  by  subtracting  the  logarithm  of  a 
baseline  image  (pre-protocol)  from  the  logarithm  of  the  remaining  images  taken  at  later  times 
while  the  gas  concentrations  were  changed  between  ambient,  hyperoxic  (carbogen  (95%  O2  and 
5%  CO2)  or  100%  O2)  and  hypercapnic  (95%  air  and  5%  CO2)  conditions.  Using  the  extinction 
coefficients  of  Hb  and  HbCL,11  and  Beer’s  Law,  we  converted  the  differential  wavelength  images 
into  changes  in  the  concentrations  of  oxygenated  and  deoxygenated  hemoglobin  (A[Hb]  and 
AfHbCL],  respectively).  “  For  these  experiments,  we  did  not  correct  for  the  scattering  in  the 
absorption  measurements,  and  thus  the  numeric  values  given  in  the  figures  are  qualitative. 

In  order  to  visualize  the  changes  in  the  images,  we  take  the  mean  and  standard  deviation 
of  each  pixel  through  time.  The  standard  deviation  image  for  [A  HbCL],  shown  in  Fig.  2, 
indicates  the  magnitude  of  the  changes  in  concentration  of  oxygenated  hemoglobin  for  the  image 
sequences,  and  demonstrates  the  remarkable  enhancement  in  differential  contrast  between 
cancerous  and  non-cancerous  tissues  in  response  to  vasoactive  gases.  As  expected,  the  majority 
of  the  response  is  contained  within  the  tumor  boundaries.  The  dynamic  response  of  the 
tumorigenic  and  non-tumorigenic  tissues  can  be  estimated  by  spatially  averaging  over  the  pixels 
outlining  the  tumor  and  the  mouse  body,  respectively,  for  each  image  throughout  the 
measurement  sequence  (Fig.  3).  The  measured  dynamics  exhibit  similar  trends  to  those 
measured  with  point-source  apparatus,  ’  namely,  that  the  concentration  of  oxy-hemoglobin 
increases  with  hyperoxia  while  the  concentration  of  deoxy-hemoglobin  decreases,  and  that  the 
total  concentration  of  hemoglobin  remains  approximately  constant.  Previous  researchers  have 
also  observed  heterogeneity  of  the  oxygenation  improvement  within  individual  tumors  by 
inserting  oxygen-measuring  probes  into  the  tumor  mass  at  different  sites.  Spatial 


3 


heterogeneities  can  be  seen  in  Fig.  2,  suggesting  that  our  non-invasive  optical  imaging  technique 
is  capable  of  probing  such  changes. 

Using  this  apparatus,  we  have  examined  many  different  tumor  models.  Each  of  the 
models  exhibits  the  same  trends  that  are  shown  in  Fig.  2  and  Fig.  3.  By  further  increasing  the 
fraction  of  CO2  delivered  to  the  mouse  to  15%,  we  have  been  able  to  see  changes  in  the  signals  at 
780  nm  and  840  nm,  which  seem  well  correlated  with  the  presence  of  O2  or  CO2,  respectively,  as 
shown  in  Fig.  4.  It  is  commonly  believed  that  vasoconstriction  results  from  breathing  pure  O2, 
whereas  a  small  fraction  of  CO2  has  a  vasodilatory  effect  on  the  mature  blood  vessels.  Recent 
work  combining  microscopy  and  MRI  have  suggested,  however,  that  both  hypercapnia  and 
hyperoxia  lead  to  vessel  constriction;  hypercapnia  leads  to  constriction  of  mature  arteries, 
whereas  hyperoxia  leads  constriction  of  all  vessels.14  We  believe  that  the  signal  changes  in  Fig.  4 
are  due  to  the  direct  vasoresponse  to  the  specific  hyperoxic  or  hypercapnic  inhalatory  conditions. 
Currently,  it  is  difficult  to  quantitatively  assess  the  values  for  the  total  concentration  of 
hemoglobin,  and  thus  it  is  difficult  to  directly  compare  our  results  with  others.  Improvements  in 
the  imaging  apparatus  in  Fig.  1  are  currently  underway,  and  should  provide  us  with  the 
We  find  that  inhalation  of  varying  levels  of  oxygen  and  carbon  dioxide  provides  strong  cancer- 
specific  contrast  in  animals.  This  approach  may  prove  useful  for  small  animal  imaging  of  cancer 
models  or  for  cancer  imaging  in  humans. 

This  research  was  supported  by  funds  from  DOD  Breast  Cancer  Research  Program,  Grant 
Number  DAMD17-02-1-0570  and  the  California  Breast  Cancer  Research  Program  of  the 
University  of  California,  Grant  Number  10EB-0168.  We  acknowledge  helpful  conversations 
with  Dr.  Mark  Dewhirst  of  Duke  University  Medical  Center. 
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Figure  Captions 


Fig.  1.  Schematic  of  the  experimental  apparatus. 

Fig.  2.  (Color  online)  Standard  deviations  of  the  individual  pixel  values  for  AfHbCF]  time-series 
of  images. 

Fig.  3.  Plot  of  A[Hb]  (a)  and  AfHbC^]  (b)  for  U87  brain  tumor  (dashed  line)  and  normal  tissue 
(solid)  in  response  to  gas  inhalation.  Gas  concentrations  were  100%  Air,  100  CF  or  95%  O2  + 
5%  CO2,  as  labeled  in  the  figure. 

Fig.  4.  Change  in  absorption  at  780  nm  (a)  and  840  nm  (b)  for  A549  lung  tumor  (dashed  line) 
and  normal  tissue  (solid)  in  response  to  gas  inhalation.  Gas  concentrations  were  100%  Air,  100 
O2  or  95%  O2  +  5%  CO2,  as  labeled  in  the  figure. 
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(54)  Title:  OPTICAL  VASCULAR  FUNCTION  IMAGING  SYSTEM  AND  METHOD  FOR  DETECTION  AND  DIAGNOSIS  OF 
CANCEROUS  TUMORS 

(57)  Abstract:  An  in  vivo  optical  imaging  system  and  method  of  identifying  unusual  vasculature  associated  with  the  angiogenic 
vasculature  in  tumors.  An  imaging  system  acquires  images  through  the  breast.  Benign,  noninvasive  oxygen  and  carbon  dioxide  are 
used  as  vasoactive  agents  and  administered  by  inhalation  to  stimulate  vascular  changes.  Images  taken  before  and  during  inhalation 
are  subtracted.  An  optical  vascular  functional  imaging  system  monitors  abnormal  vasculature  through  optical  measurements  on  oxy- 
and  deoxyhemoglobin  during  inhalation  of  varying  levels  of  02  and  C02.  The  increase  in  contrast  between  tumor  (cancerous)  and 
normal  (noncancerous)  tissue  is  dramatic,  facilitating  accurate  early  detection  of  cancerous  tumors  and  improving  sensitivity  and 
specificity  (lower  false  negative  and  false  positive  rates).  The  invention  is  useful  in  mammography,  dermatology,  prostate  imaging 
and  other  optically  accessible  areas. 


